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Abstract 



A new package, DISPred, is described. The package can be used to calculate e ± p deep 
inelastic scattering cross sections at Born level in Electroweak theory and at both leading 
and next-to-leading order in QCD. 

> 

o\ : 1 Introduction 

lyt . The package DISPred arose as a result of the need to produce predictions of deep inelastic 

scattering (DIS) electron-proton cross sections at next-to-leading order (NLO) in QCD 
for comparison to data in ZEUS publications. In the current version (1.0) predictions for 
the following are available at both leading order (LO) and NLO in QCD: 

• the reduced/ double-diff cross sections for neutral current (NC) and charged current 
(CC) DIS 

X ■ 

$—i ' • the total cross section a tot for NC and CC DIS 

• differential cross sections M^ g, g for NC and CC DIS 

for collisions of unpolarized beams of electrons (e ± ) and protons. 

DISPred has been tested with LHAPDF 5.8.2 [lj and can produce predictions for the 
ZEUS- JETS [2\ or HERA0.1 parton distribution functions (PDFs) in LHAPDF in the 
.LHpdf format, or any other PDF within LHAPDF in the .LHgrid format. 

It produces output in ascii text format and can also produce histograms and graphs 
in ROOT-based [3] formats. 

2 Leading Order Calculation 

The LO QCD, Born- level electroweak, cross section is calculated according to the formu- 
lation given by Devenish and Cooper-Sarkar[4J. 
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2.1 Reduced and Double Differential Cross Sections 
2.1.1 NC DIS 

The double differential cross section in NC scattering is: 



.,=„ 



W = 2 ^ [F ^ Q2) - T + F ™^ * F^ N0 (*,Q 2 )]. (i) 

Where, as is conventional, Q 2 is the virtuality of the exchanged boson, x is the mo- 
mentum fraction of the struck parton in the infinite proton- momentum, F^ are structure 
functions defined later, Y± = 1 ± (1 — y) 2 and y is the inelasticity of the electron. 

For the leading order calculation the structure functions are defined as follows: 

if C =Y,A°(Q 2 )(xq l (x,Q 2 )+xq l (x,Q 2 )); (2) 

i 

F L NC = 0; (3) 

^3 NC = E B°(Q 2 )(x qi (x, Q 2 ) - xq t (x, Q 2 )); ; (4) 

i 

where Ai and Bi can be expressed in terms of the NC vector and axial-vector elec- 
troweak couplings to the quarks (electron) v i (v e ) and a^ (a e ) and quark charge e^ as 



and 



and 



4° = e 2 - 2e l v i v e P z (Q 2 ) + (v 2 + a 2 e )(v 2 + a 2 )P 2 z (Q 2 ) (5) 

£° = -2e iai a e P z (Q 2 ) + 4a iVl v e a e P 2 (Q 2 ) (6) 



p ^ 2 ) = ^fu|(^W) < 7 » 



The reduced cross section for NC scattering is: 






*ncT = l^ C (x,Q 2 ) " ^F^(x,Q 2 ) T ^xF 3 NC (x,g 2 )]. (8) 



2.1.2 CC DIS 

The double differential cross section in CC scattering is: 
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Where Myy is the mass of the W boson, Gp the Fermi coupling constant and the Fp c 
are defined at LO in QCD a^]: 

(10) 
(11) 
(12) 
(13) 

which leads to the following expressions for the reduced cross sections: 



of+ c = x[(l-y) 2 (d + s)+u + c] (14) 

a c e ? = x[(l-y) 2 (d + s)+u + c] (15) 

2.2 Single Differential and Total Cross Sections 

In order to calculate the single differential and total cross sections, the expression for the 
double differential cross sections is integrated over the allowed regions of Q 2 , x and y using 
the VEGAS [5] algorithm as implemented in the GNU Scientific Library [6j. The number 
of calls used in VEGAS may be specified via the control cards. Differential cross sections 
may also be calculated at a "point"; in this case the width of the bin which contains 
the point is multiplied by a predetermined factor (which can be chosen in the control 
cards) to provide an approximate calculation. If no input points are specifed via cards, 
the option "AUTO" may be chosen which makes the program DISPrediction calculate 
the differential cross sections at the centroid of the bin. 



3 Next-to-Leading Order Calculation 

3.1 Reduced and Double Differential Cross Sections 

The implementation of Qcdnum 16.13 |7J included in LHAPDF is used in DISPred to 
evaluate structure functions i*2, Fl and F% at NLO in QCD. The prescription used by 
the ZEUS collaboration for the ZEUS JETS fit [SI [2] has been adopted. As such DISPred 
can perform the QCD evolution for the ZEUS-JETS and ZEUS-S fits and use the . LHpdf 
format files from LHAPDF for this. In the case of other PDFs DISPred can fill a Q 2 ,x 
grid for Qcdnum using the values from the tt .LHGrid file. The structure functions are 
then generated from this grid. All other aspects of the reduced cross section cross section 
are the same as for the leading order case. 

Predicted NLO reduced cross sections in CC DIS made using DISPred for e + p col- 
lisons with proton energy 920 GeV and positron energy 27.56 GeV are shown in figure [TJ 
Predictions are shown for the PDF sets ZEUS- JETS 0, MSTW08 [9], CTEQ66 [TO] and 
HERAPDF1.0 [11 j . In addition the uncertainties for the ZEUS-JETS predictions are 
shown as a yellow band. 



^Jn the expressions shown above the small, Cabbibo-suppressed, contribution from the 6-quark is neglected, 
it is however included in the calculation made by DISPred. In vl.O DISPred is only suitable for use for the 
HERA energy regime and so top quark contributions are not included 
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Figure 1: Predictions of a at NLO in QCD for CC DIS in e + p collisions with proton energy 
920 GeV and positron energy 27.56 GeV. 



3.2 Single Differential and Total Cross Sections 

The single and total cross sections at NLO are calculated in precisely the same manner as 
for the LO calculations except that the NLO expressions for the structure functions are 
used. 



4 Installation and Usage 

4.1 Requirements 

The code has been tested on GNU/Linux systems and as such the following packages are 
required for installation: 

• GNU autoconf automake 

• The GNU scientific library 

• LHAPDF 

• ROOT 

4.2 Installation 



Tarballs of the package may be downloaded from http : //www . hepf orge . org/downloads/dispred. 

After downloadig the tarball, the package may be installed with: 

tar -zxvf DISPred-l.O.tgz 

cd DISPred 

./configure — pref ix=<installation directory> 

make 

make install 

This assumes that root-conf ig and lhapdf-conf ig are already in your path. Do not 

choose <installation directory> to be the same as the location of the expanded tar 

file. It is now possible to run the example program or to include the DISPred libraries in 

your own code. 

4.3 Example Program: DISPrediction 

Within <installation directory>/bin there is a program DISPrediction which can 
be used to produce predictions for ep DIS cross sections. This program takes as input a 
cards file e.g.: 
DISPrediction example . cards 

An example set for testing is available in the example subdirectory of the tarball. The 
available options for the cards file are summarised in table [TJ 

4.4 DISPred Library and Classes 

The DISPred packagedprovides a library as well as the DISPrediction executable. This 
library makes it easy to construct programmes that calculate DIS cross sections. An 
example of a simple programme is in fact DISPrediction itself, which is very short: 



Type 


Name 


Default Value 


Meaning 


int 


ElectronCharge 


1 


choose e + or e~ 


int 


VegasCalls" 


50000 


Calls to Vegas for integration 


string 


PDFSetFileName 


"ZEUS2005_ZJ.LHpdf" 


Name of PDF file to use 


int 


PDFSubSet 





Subset of PDF to use 


string 


QCDCalculationLevel 


"LO" 


Can be "LO" or "NLO" 


string 


DISProcess 


"NC" 


Choose CC or NC 


double 


ZBosonMass 


91.1876 


M z (GeV) 


double 


WBosonMass 


80.398 


M w (GeV) 


double 


AlphaEM 


7.297352570x10-3 


«EM 


double 


GFermi 


1.1663710- 5 


G F 


double 


TopMass 


171.2 


M t 


double 


BottomMass 


4.20 


M b 


double 


Vub 


41.2e-3 


CKM V ub 


double 


Vcb 


3.93e-3 


CKM V cb 


double 


Sin2ThetaW 


0.22308 


sin 2 9 W 


double 


Sin2ThetaC 


0.05 


sin 2 9 C 


double 


CouplingVu 


0.203 


v u SM=0.5-4*sin 2 6W 


double 


CouplingVd 


-0.351 


v d SM= -0.5 + 2 * sin 2 9 W 


double 


CouplingVe 


-0.00538 


v e SM=-0.5 + 2 * sin 2 9 W 


double 


CouplingAu 


0.5 


a u 


double 


CouplingAd 


-0.5 


ad 


double 


CouplingAe 


-0.5 


a e 


string 


ReducedCrossSection 


"OFF" 


can be OFF or ON 


string 


ReducedCrossSectionBins 


"q2xpoints.dat" 


file containing points for a 


double 


DiffBinPointScale 


le-6 


Fraction of bin width for -M^ etc. 


string 


DSigmaDQ2 


"OFF" 


can be OFF or ON 


string 


DSigmaDQ2Bins 


"q2bins.dat" 


file containing bins for 


string 


DSigmaDQ2Points 


"AUTO" 


file with points for J% 


string 


DSigmaDX 


"OFF" 


can be OFF or ON 


string 


DSigmaDXBins 


"xbins.dat" 


file containing bins for ^ 


string 


DSigmaDXPoints 


"AUTO" 


file with points for 


string 


DSigmaDY 


"OFF" 


can be OFF or ON 


string 


DSigmaDYBins 


"ybins.dat" 


file containing points for ^- 


string 


DSigmaDYPoints 


"AUTO" 


file with points for 


double 


Q2Min 


0.0 


minimum Q 2 


double 


Q2Max 


100000.0 


maximum Q 2 


double 


XMin 


0.0 


minimum x 


double 


XMax 


1.0 


maximum x 


double 


YMin 


0.0 


minimum y 


double 


YMax 


1.0 


maximum y 


double 


ELepton 


27.5 


Electron beam energy 


double 


EProton 


920.0 


Proton beam energy 


string 


ROOTOutputFile 


"DISPredOut.root" 


Output file for ROOT 



Table 1: Available control cards for DISPrediction. 
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#include <iostream> 
#include "DISPredictor .h" 

using namespace DISPred; 

int main (int argc, char **argv){ 

std::cout « "DISPrediction vl.O - 31 Mar 2010" « std::endl; 

// Create instance of DISPredictor 

DISPred: : DISPredictor *DISPred= DISPred: : DISPredictor :: Instance () ; // initialise 

from control cards provided via command line 

DISPred->Initialise(argc,argv) ; DISPred->CalculateCrossSections() ; 

DISPred->PrintResults() ; 

std::cout << "DISPrediction vl.O - Run finished Succesfully" << std::endl; 

DISPred->WriteOutput() ; 

return 0; 

} 

All classes are part of the name space DISPred. 

4.4.1 The ControlCards Class 

The control cards class is used to handle configuration options that can be read in from 
a text file. It is implemented as a singleton class. Available methods for the class are 
detailed below. 

ControlCards* Instance (): Returns a pointer to the instance of control cards. 

void AddCardDouble (const std::string key, const double defval): Definesacard 
with name key and with a default double precision value defval. 

void AddCardlnt (const std: : string key, const int defval): Defines a card with 
name key and with a default value defval which is an integer. 

void AddCardString (const std::string key, std::string defval) : Definesacard 
with name key and with a default value defval which is a string. 

void AddCardVector (const std::string key, const std: :vector<double> defval): 
Defines a card with name key and with a default value defval which is a vector of 
double precision values. 

int readKeys (const char* fileName): Reads in card values from the file with name 
f ileName. 

double fetchValueDouble (const std: :string& key): fetch the value of card key. 

int f etchValuelnt (const std::string& key): fetch the values of card key. 

std: : string f etchValueString(const std: :string& key): fetch the value of card 
key. 

std: :vector<double> fetchValueVector (const std::string& key): fetch the val- 
ues of card key. 

void printCardsQ : Print current card values to stdout. 



4.4.2 The DISPredictor class 

The DISPredictor class is a singleton class that is the workhorse of DISPred. It has 
many public methods. 

static DISPredictor* Instance (): returns the instance of DISPredictor. 

void Initialise (int my_argc, char **my_argv): Initialise DISPredictor based on a 
cards file name which can come directly from stdin. 

void CalculateCrossSectionsO: Calculate cross sections as configured in the cards. 

void InitPDF(int subset): Initialise the chosen PDF set. 

void PrintResults(): Print results to stdout. 

void WriteOutput (): Write the output rootfile. 

double CalculateReducedCrossSection(double x, double q2): Calculate a NC DIS 
reduced cross section. 

double CalculateCCReducedCrossSection(double x, double q2): Calculate a CC 
DIS reduced cross section. 

double CalculatePropagator (double q2, double x): Calculate the NC propagator. 

double CalculateCCPropagator (double q2, double x): Calculate the CC propaga- 
tor. 

double CalculateDSigmaDQ2 (double q2min, double q2max): Calculate -77%. 
double CalculateDSigmaDX (double xmin, double xmax): Calculate 4^. 



double CalculateQ2DSigmaDQ2 (double q2min, double q2max): Calculate Q 2 d ' 7 



double CalculateXDSigmaDX (double xmin, double xmax) :Calculate x4^z- 

double CalculateDSigmaDY (double ymin, double ymax): Calculate -£-. 

double CalculateYDSigmaDY (double ymin, double ymax): Calculate y^-- 

double S(): Return the centre-of-mass energy squared. 

4.4.3 The RedSigmaPoint class 

The RedSigmaPoint class is a simple class for storing information about double-differential 
cross sections points. For each point the Q 2 (-q2), 2;(_x), a(_redsigma) and dQ i d (_d2sdq2dx) 

RedSigmaPoint (double q2, double x): constructor that creates a point with _q2=q2 
and _x=x and other values 0. 

RedSigmaPoint (double q2, double x, double redsigma): constructor that creates a 
point with _q2=q2 and _x=x, _redsigma=redsigma and _d2sdq2dx=0. 

RedSigmaPoint (double q2, double x, double redsigma, double d2sdq2dx ): con- 
structor that creates a point with _q2=q2 and _x=x, _redsigma=redsigma and 
_d2sdq2dx=d2sdq2dx. 

RedSigmaPoint (): Constructor with all vlaues set to 0; 

double Q2(): returns _q2. 



double X(): returns _x. 

double RedSigma(): returns _redsigma. 

double D2sDQ2Dx(): returns _d2sdq2dx. 

void SetRedSigma(double reduced): Set _redsigma. 

void SetD2sDQ2Dx (double reduced): Set _d2sdq2dx. 

void Print (): Print out information. 

void PrintShort (): Briefly print out information. 

4.4.4 The RedSigmaGrid class 

The RedSigmaGrid class inherits from a std: : vector<RedSigmaPoint>. With the fol- 
lowing extra methods: 

void Print (): Print out information. 

void PrintShort (): Briefly print out information. 

4.4.5 The Dif f SigmaPoint class 

The Dif f SigmaPoint class is a simple class for storing information about single-differential 
cross sections at a point. A point in the variable of choice called _var and the differential 
cross section _dif f sigma are stored. The following public methods are available. 

Dif f SigmaPoint (): Default constructor, sets _var to 1.5 and _diff sigma to 0. 

Dif f SigmaPoint (double var): Constructor that creates a point at _var= var with 
_diff sigma=0. 

Dif f SigmaPoint (double var, double dif f sigma): Constructor that creates a point 
at _var= var with _diffsigma=diff sigma. 

void Print (): Print out information. 

void PrintShort (): Briefly print out information. 

double Var(): returns _var. 

double DiffSigma(): returns _diff sigma. 

void SetDiffSigma(double dif f sigma): sets _diff sigma to dif f sigma. 



5 Root Output 



When DISPred produces an output root file, then a TTree and several histograms and 
graphs are produced. 

5.1 Root TTree 

A TTree called ReducedCrossSections is produced. The variables in this tree are listed 
in table [2j 



Type 



Variable Name 



Description 



An integer giving the ID of the point 

The Q 2 of the point 

The x of the point 

The double-differential cross section 
The reduced cross section 
The d PDF at this point 
The d PDF at this point 
The u PDF at this point 
The u PDF at this point 
The s PDF at this point 
The s PDF at this point 
The c PDF at this point 
The c PDF at this point 
The b PDF at this point 
The b PDF at this point 



int 
double 
double 
double 
double 
double 
double 
double 
double 
double 
double 
double 
double 
double 
double 



point 

Q2 

x 

ddiff sigma 

redsigma 

d 
dbar 

u 
ubar 

s 
sbar 

c 
cbar 

b 
bbar 



Table 2: Tree variables in the root output file. 



5.2 Root Histograms 

Six TH1D objects are produced: 
DSigmaDQ2: Binwise -M^; 
DSigmaDX: Binwise *jf ; 
DSigmaDY: Binwise <j-. 
Q2DSigmaDQ2: Binwise Q 2 J%; 
XDSigmaDX: Binwise x^; 
YDSigmaDY: Binwise y^-. 

5.3 Root Graphs 

Three TGraphAsymErrors objecs are produced: 
GraphDSigmaDQ2: Pointwise -M?; 

GraphDSigmaDX: Pointwise ||; 
GraphDSigmaDY: Pointwise 4f. 



6 Summary 



This manual for the DISPred package vl.O has outlined the features currently imple- 
mented together with a simple example programme that will make predictions for DIS 
cross sections in ep scattering. The code and most up-to-date information are hosted by 



hepforge at: http://projects.hepforge.org/dispred/ 



10 



Acknowledgements 



The author wishes to thank C. Gwenlan for help with QCDNUM and cross checks of 
the results from DISPred, M. Sutton for testing the code, A. Tapper for providing code 
that is used for the implementation of control cards and R. Ciesielski, A. Cooper-Sarkar, 
K. Oliver, E. Tassi and M. Turcato for feedback on the results. 

References 



[i] 

[2] 

[3] 

[4] 

[5] 
[6] 

[7] 
[8] 

[9] 
[10] 
[111 



M. R. Whalley, D. Bourilkov, and R. C. Group, The Les Houches Accord PDFs 



(LHAPDF) and Lhaglue, hep-ph/0508110. 



ZEUS Collaboration, S. Chekanov et al., An NLO QCD analysis of inclusive 
cross-section and jet- production data from the zeus experiment, Eur. Phys. J. C42 
(2005) 1-16, |hep-ph/0503274| . 

R. Brun and F. Rademakers, ROOT: An object oriented data analysis framework, 
Nucl. Instrum. Meth. A389 (1997) 81-86. 

R. Devenish and A. Cooper-Sarkar, Deep inelastic scattering, . Oxford, UK: Univ. 
Pr. (2004) 403 p. 

G. P. Lepage, A New Algorithm for Adaptive Multidimensional Integration, J. 
Comput. Phys. 27 (1978) 192. 

M. Galassi et al, GNU Scientific Library Reference Manual, . (3 Ed.) ISBN 
0954612078. 

M. Botje, A QCD analysis of HERA and fixed target structure function data, Eur. 
Phys. J. C14 (2000) 285-297, |hep-ph/99 12439] . 

ZEUS Collaboration, S. Chekanov et al, A ZEUS next-to-leading- order QCD 
analysis of data on deep inelastic scattering, Phys. Rev. D67 (2003) 012007, 
|hep-ex/0208023]. 

A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, Parton distributions for 
the LHC, iarXiv: 0901. 00021 

W. K. Tung et al., Heavy quark mass effects in deep inelastic scattering and global 
QCD analysis, JHEP 02 (2007) 053, |hep -ph/06lT254] . 

HI and ZEUS Collaboration, F. D. Aaron et al., Combined Measurement and 
QCD Analysis of the Inclusive ep Scattering Cross Sections at HERA, JHEP 01 
(2010) 109, [arX iv: 0911 .0884) . 



11 



